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i!B ■ ABSTRACT 

^ I Ring waves of star formation, propagating outwardly in the galactic disks, leave 

■ chemical abundance gradients in their wakes. We show that the relative [Fe/0] abun- 

^ ■ dance gradients in ring galaxies can be used as a tool for determining the role of the SNIa 

explosions in their chemical enrichment. We consider two mechanisms which can create 
outwardly propagating star forming rings in a purely gaseous disk — a self-induced wave 
^ ' and a density wave, and demonstrate that the radial distribution of the relative [Fe/0] 

. abundance gradients does not depend on the particular mechanism of the wave forma- 

tion or on the parameters of the star- forming process. We show that the [Fe/0] profile 
is determined by the velocity of the wave, initial mass function, and the initial chemical 
composition of the star-forming gas. If the role of SNIa explosions is negligible in the 
chemical enrichment, the ratio [Fe/0] remains constant throughout the galactic disk 
with a steep gradient at the wave front. If SNIa stars are important in the production 
of cosmic iron, the [Fe/0] ratio has gradient in the wake of the star-forming wave with 
the value depending on the frequency of SNIa explosions. 

Subject headings: GALAXIES: structure — ISM: kinematics and dynamics, structure 
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1. INTRODUCTION 

Theoretical modeling of chemical abundance gradi- 
ents in galaxies provides an important tool for under- 
standing galactic evolution. However, chemical evolu- 
tionary models applied for the dynamics of galaxies on 
a Hubble time scale encounter some principal difficul- 
ties (Shore & Ferrini 1995). The chemical evolution of 
galaxies is usually considered in the one-zone approx- 
imation with a postulated accretion rate and some 
assumed law of star formation. These models do not 
take into account global hydrodynamic processes of 
the redistribution of the matter in the galactic disks. 
The spiral density waves, viscosity, and the disk-halo- 
disk circulation work in different ways to redistribute 
matter within the galactic disks during the galactic 
evolution, which considerably complicates the prob- 
lem. 

We discuss one class of galaxies where the modern 
theories of the chemical evolution can be successfully 
applied. These are the starburst galaxies with the 
large scale rings of an enhanced star formation. These 
galaxies are believed to be the result of recent galactic 
collisions, and result from the passage of a compact 
companion galaxy through the disk of a galaxy along 
its rotation axis (Lynds & Toomre 1976). Such col- 
lisions produce the wave of enhanced star formation, 
and are responsible for the nature of the large scale 
rings of the young massive stars observed in some 
galaxies. 

The time-scale of bursts of star formation in the 
galaxies similar to the Cartwheel is much shorter com- 
pared to the Hubble time, making secular hydrody- 
namic processes unimportant. The chemical evolu- 
tion in the fast waves of star formation is governed 
therefore by a few parameters which can be measured, 
or at least estimated, and ring galaxies provide a pos- 
sibility to determine the relative role of SNIa and SNII 
in the enrichment of the interstellar matter by heavy 
elements. 

It is known that there are two physically different 
sources of heavy element production. The a-elements 
and oxygen arc mainly produced in massive stars ex- 
periencing SNII explosions (Woosley & Weaver 1995), 
and iron peak elements are substantially contributed 
by SNIa which are assumed to be the end product of 
close binary evolution (Tsujimoto et al. 1995). The 
relative role of these mechanisms is somewhat unclear 
due to uncertainties in knowledge of the input param- 
eters, and uncertainties of stellar evolutionary models 



(Ishimaru & Arimoto 1997; Gibson, Loewenstein & 
Mushotsky 1998). 

The clock of the SNIa explosions and thus the en- 
richment due to SNIa is substantially delayed with 
respect to SNII. Such a delay in iron chemical enrich- 
ment compared to the a-elements will reveal itself in 
the star-forming ring galaxies as a relative abundance 
[Fe/a] radial gradient. The value of this gradient de- 
pends only on a few parameters such as the velocity 
of the wave, the slope of the initial mass function, and 
on the fraction of close binary systems producing iron 
via type SNIa explosions. The number of "free" pa- 
rameters is hence drastically reduced, and ring galax- 
ies give us unique opportunity to make observational 
conclusions about the role of SNIa explosions in chem- 
ical evolution. 

2. CHEMICAL EVOLUTION MODEL 

As in our previous paper (Korchagin et al. 1998), 
we discuss two possible mechanisms of the forma- 
tion of large-scale star-forming rings. The conven- 
tional scenario assumes that rings originate in the 
direct collisions of disk galaxies with their compan- 
ions. Instead, we consider the possibility that galactic 
collision only plays the role of a "detonator" stimu- 
lating further self-propagating star formation. This 
model assumes that the wave is a self-organized phe- 
nomenon, and is similar to the "fire in the forest" 
models discussed earlier. 

The most observationally studied ring galaxy is the 
Cartwheel, and we adopt therefore parameters of the 
model obtained for this archetype galaxy. The ve- 
locity of the wave was chosen equal to 90 km/s - 
the value, which gives the best fit to the optical sur- 
face brightness gradients observed in the Cartwheel 
galaxy (Korchagin et al. 1998). In both models we 
use the Salpeter IMF (a = 1.35) with masses of stars 
located in the interval 0.5 < TOs < 50 M©. It is 
remarkable however, that the ratio of iron to oxygen 
production in the wave does not depend on the pa- 
rameters determining the rate or the efficiency of star 
formation. This ratio is mainly determined by the 
IMF slope, and by the heavy element output of the 
SN explosions. This fact makes the theoretical pre- 
dictions of the radial dependence [Fe/0] ratio quite 
robust. 

The mathematical formulation of both models, as 
well as the parameters prescribing the rate of star 
formation can be found in Korchagin et al. (1998). 
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The basic ideas of the chemical evohition model 
which we incorporate here to describe abundance gra- 
dients in the starburst ring galaxies were outlined by 
Tcnorio-Taglc (1996). The process of clement enrich- 
ment in galaxies can be summarized as follows. The 
freshly formed elements ejected by the SN explosions 
fill the hot interiors of the bubbles, which were cre- 
ated earlier by the stellar winds, and/or by the re- 
lease of energy in a solitary, or multiple SN explosions. 
The new elements arc mixed with pre-existing matter 
through evaporation of a part of the cold shell sur- 
rounding the bubble, before they are expelled into the 
halo of the galaxy in a hot, SN-driven winds ("chim- 
neys"). Hot enriched gas eventually cools and "rains" 
back on the galactic disk completing the circle of the 
element enrichment in a new event of star formation. 

Recent observations support this picture. Kobul- 
nicky (1999) did not find any evidence for metal abun- 
dance enhancement in the vicinity of young star clus- 
ters, concluding that the freshly ejected materials are 
stored in a hot, 10® gas phase. Elmegreen (1997), 
considering processes of mixing and contamination of 
interstellar gas after passage of a spiral density wave, 
came to a similar conclusion. He finds that there is 
no element mixing on time scales of a few x 10^ years. 

The following processes determine spatial distribu- 
tion of heavy elements in the wake of a star-forming 
wave. 

a) Stellar Evolution: The new-born stars "inherit" 
heavy element abundances of the pre-existing gaseous 
disk, and hence their atmospheres are not expected to 
contain the heavy elements produced in the current 
starburst. Gas released by the low-mass single stars 
during their red giant phase and the SN explosions 
with mass m < 11 is assumed to have oxygen and 
iron abundances equal to that at the birth epoch of 
the stars. We take this abundance as to be one fifth 
of the solar. The fraction of mass returned to the 
interstellar medium by stars with to < 11 during their 
evolution is taken from Koppen & Arimoto (1991). 

Stars with masses m > 11 Mq experience SNII ex- 
plosions and release freshly formed elements. We use 
in our simulations the oxygen and iron yields taken 
from models of Woosley & Weaver (1995) and Nomoto 
et al. (1997). 

The SNla result from C-dcflagration of white dwarfs 
in binary systems with masses of progenitors < 8Mq. 
The nature of the SNI progenitors is rather unclear, 
and different scenarios give the lifetimes of SNIa pro- 



genitors in the interval ^ 0.1 — 0.6 Gyr (Branch 1998). 
In our simulations, we consider therefore three differ- 
ent mass intervals for the SNIa progenitors. Namely, 
we assume that SNIa are produced in binaries with 
primary stars of 2.5-8 Mq, 2.5-6 Mq, and 2.5-3.5 
Mq. Hence the first supernovae explode at 4 x 10^, 
7x 10^ and 2.66 x 10^ years after the start of the burst, 
for the three mass intervals. In the latter case, SNIa 
become unimportant in the heavy element enrichment 
by fast ring waves of star formation such as inferred 
in Cartwheel. We adopt the nucleosynthesis output 
for SNIa from the updated W7-model of Thielemann, 
Nomoto & Hashimoto (1993) with fraction of bina- 
ries fsNi equal to 0.04 — the value which Kobayashi 
et al. (1998) found to be adjusted to reproduce the 
chemical evolution in the solar neighborhood. 

b) Evaporation: There are two sources of mixing 
of the chemically unprocessed gas with the freshly 
formed elements inside the hot bubble. The first 
mechanism is the classical thermal conductivity, which 
results typically a few thousand solar masses being 
evaporated from the shell surrounding the hot bub- 
ble (e.g. Shull & Saken 1995). A larger potential 
source of mixing, but far more uncertain, is the photo- 
evaporation or shock ablation of the ambient clumpy 
medium. We consider therefore the following models 
of the chemical enrichment in the Cartwheel. In the 
first model, we assume that the freshly formed ele- 
ments are instantaneously mixed with the remaining 
ambient gas, or equally, that all remaining ambient 
gas is transformed into a hot phase by the star for- 
mation process. In the second model, we consider 
the conductive evaporation of shells to be the source 
of mass in the bubble. For comparison, we also com- 
puted the [Fe/O] gradients for pure SN ejecta without 
mixing. 

Using expression for the conductive mass-loss of 
the shell in the adiabatic bubble (Shull & Saken 1995) 
and assuming 10^^ erg of energy released by each SN 
and the shell expansion time of 10^ yrs, one can write 
the expression for the SN-evaporated mass in units 
lO^M©, 10® yr and 1 kpc as: 

Me^ap = (0.7 X 10-4)no(r, t)-^/^^K^/' (1) 

Here, no{r,t) is the density of the ambient gas, and 
K < 1 is a scaling factor. 

c) Cooling: The hot processed gas streams out in 
the halo, and the subsequent abundance gradients 
are determined by the concurrent processes of cool- 



3 



ing and diffusion. Using Raymond's ct al. (1976) 
cooling coefficient we find that the coohng time of 
hot X-ray emitting gas with the temperature ~ few 



xlO^K and tlio density 10" 



10- 



is about 



a few x(10^ — 10^) yrs which is large, or comparable 
to the dynamical time of the star burst. At the tem- 
perature interval typical for the X-ray emitting halo 
gas, the heavy elements, particularly iron, dominate 
the cooling process, accounting from half to three- 
forth of the total cooling rate (Raymond ct al 1976). 
The heavy element abundances in the Cartwheel are 
deficient by factor of ten as compared to the Orion 
Nebula in our galaxy (Fosbury & Hawarden 1977), 
which will increase the estimated cooling time. We 
assumed therefore that cooling is unimportant dur- 
ing the starburst time in the Cartwheel. 

d) Diffusion: The heavy ion admixtures diffuse 
into a hot hydrogen plasma of halo, decreasing pos- 
sible abundance gradients produced by the SN explo- 
sions. The diffusion of light ions is effective on the 
time scale of the starburst (Tenorio-Tagle 1996), but 
is ineffective for heavy, highly ionized ions. In plasma 
with temperature of a few xlO^ K, the oxygen will 
be ionized up to OVI, and iron up to FeXI-FelV. The 
mean free paths of oxygen and iron in plasma with 
temperature 5 x 10^ K and density 10"'^ will be about 
50pc and 12pc correspondingly, which gives values of 
the dimensionless diffusion coefficients Dq « 5 x 10~^ 
and Dpe « 6 x 10~*. Our simulations show, that with 
such values, the diffusion of heavy elements is unim- 
portant on the time-scales of passage wave front in 
the disk. 

e) Grain Formation: For the pressure typical in a 
SN-processed gas, the grain formation is effective for 
temperatures below 1000 K (Gail & Sedlmayr 1985). 
Hot halo gas will preserve therefore the abundance 
compositions which were formed during the ejecta. 
However, if the hot gas had enough time to cool, iron 
will be selectively depleted with respect to oxygen due 
to grain formation. 

Assuming that the SN ejecta mix with all the re- 
maining ambient gas, we can write the equation de- 
scribing the spatial distribution of the heavy elements 
in the wake of the star-forming wave: 

'""^^^'^^^ -Z05(M) + A^ 



dt 



,)(l)(rn)R,{vi)dm (2) 



Here B{r, t) denotes the birth rate of stars as a func- 
tion of time and radius of the galactic disk, Mj is the 



surface density of an i-th element, = Mf /Mqas is 
the initial abundance of the i-th element in the star- 
forming gas (/>(m) is the initial mass function with 
minimum of maximum masses Af„ji„ and M^^ax re- 
spectively, is the diffusion coefficient of i-th ele- 
ment, and Mgas is the initial constant surface density 
of gas. The coefficient Ri{m) gives ratio of mass of 
the ejected element to the mass of a star to. 

In a similar way, we can write the equation gov- 
erning the chemical evolution behind the wave when 
the stellar ejecta are mixed with gas thermally evap- 
orated from the shells with its mass prescribed by the 
equation (1). 

Once the IMF is fixed, the birth and death rates 
can be determined at each radial grid zone making 
use of the lifetimes of stars. Then the heavy element 
distributions can be computed by means of equations 
(1) and (2). 

3. Results 

Figure 1 illustrates the [Fe/O] enrichment pro- 
duced by the wave of star formation for the three 
different scenarios of the element mixing and three 
different mass intervals of the pre-SNIa stars. The 
position of the wave taken at time t = 2AQMyr cor- 
responding to the present location of the outer ring 
in the Cartwheel galaxy. The left frames show the 
abundance gradients for the Nomoto et al. (1997) 
SNII outputs, and the right frames present the re- 
sults of simulations for the Woosley & Weaver (1995) 
SNII iron and oxygen outputs. 

The bottom frames of Figure 1 show the [Fe/O] 
profiles if masses of pre-SNIa stars 2.5-3.5 Mq. In 
this case, all SNIa progenitors have lifetimes longer 
than the time of the wave propagation , and the role 
of SNIa stars in the enrichmcnit is negligible. The solid 
lines of the Figure 1 show the [Fe/0\ radial depen- 
dence when hot SN ejecta are mixed with the part of 
ambient gas evaporated from the SN shells. For com- 
parison, we also plot the abundance gradients when 
the freshly synthesized elements are instantaneously 
mixed with the all ambient gas (dashed lines) , and the 
abundance gradients produced by the pure SN ejecta 
(dotted lines). 

Figure 1 depicts more sharp growth of the oxygen 
abundance as compared to the iron one in front of the 
wave of star formation. Such behavior directly follows 
from the yields of iron and oxygen by the SNII explo- 
sions. In the model of Woosley & Weaver (1995) the 



4 



oxygxni yield increases with stellar mass until the pro- 
genitor mass reaches approximately 40 Mq, and then 
tends to saturate. Supernova models of Nomoto et 
al. (1997) give rapid increase of oxygen yield with the 
mass of progenitor. The iron yield on the contrary de- 
creases with stellar mass in both the models. There- 
fore, SNII with moderate masses of progenitor are 
responsible mainly for the iron production, whereas 
SNII with progenitors of higher mass produce mainly 
oxygen. 

Large [Fe/O] gradients in the region of the wave 
of star formation are changed to slower growth of the 
[Fe/O] ratio towards the center of galaxy when metal 
enrichment due to SNIa events start to become im- 
portant. The top and intermediate frames of Figure 
1 show the radial profiles of [Fe/O] for the SNIa pro- 
genitors with masses 2.5-8 Mq and 2.5-6 Mq, respec- 
tively. The "pollution" by the SNIa hot ejecta causes 
the increase of the [Fe/O] gradient towards the center 
of the disk. The effect is obviously most noticeable 
for the pure SN ejecta (dotted lines). The element 
mixing diminishes the effect, which however remains 
about 0.05 dex in the Fe/O ratio in both SNII output 
models and both models of mixing of newly-formed 
elements. 

The radial abundance profiles of oxygen and iron 
obtained for the density wave are much similar to 
those obtained for the induced wave of star forma- 
tion. The radial distributions of relative abundance 
[Fe/O] are qualitatively independent of the particular 
model of star formation and the uncertain parameters 
such as the rate of star formation and the "efficiency" 
of star formation. This makes our conclusions robust, 
and the relative abundance gradients [Fe/a] might be 
used as a tool for determining the relative role of SNIa 
in the chemical enrichment. 

4. DISCUSSION 

The nebular abundances measured by Fosbury & 
Hawarden (1977) in the Cartwheel pertain to the 
heavy element abundances before the passage of the 
wave. Therefore, observations in the inner ring re- 
gions are required to test our predictions. 

Interstellar abundances can be inferred however us- 
ing X-ray emission from hot halo gas. The launch of 
the Advanced X-ray Astrophysics Facility (AXAF) 
with its spatial resolution of ~ 0."5 would allow 
to measure the heavy element gradients in the ring 
galaxies of star formation. 



If the hot gas has enough time to cool after the 
passage of the wave of star formation, and the newly- 
formed elements reside in the warm ionized or the cold 
phases of interstellar medium, then abundances can 
be inferred using metallic absorption lines towards 
background objects. The absorption lines of most 
atoms found in the interstellar medium occur at ultra- 
violet wavelengths, and hence abundance determina- 
tions require satellite observations. With the launch 
of Hubble Space Telescope interstellar abundances are 
now available for a number of atoms, including iron 
and the a elements O, N, C, Mg etc. (Savage & Sem- 
bach 1996). 

In external galaxies, the measurement of interstel- 
lar abundances depends on the availability of back- 
ground sources. The most widely used background 
sources are quasars, using which abundances in halos 
of galaxies and in the intervening clouds have been 
determined (e.g. Morton ct al. 1980). In a recent 
study, Gonzalez-Delgado et al. (1998) have detected 
interstellar absorption lines towards the bright star- 
burst nucleus of NGC 7714 in the ultraviolet spectrum 
taken with the Hubble Space Telescope. Hence, when 
background sources are available the present technol- 
ogy allows for the measurement of interstellar absorp- 
tion lines. In ring galaxies, the measurement of inter- 
stellar abundance thus depends on the inclination of 
the ring to the line of sight. In cases, where the ring is 
favorably aligned, the spectrum of center regions and 
the ring HII regions are likely to contain interstellar 
absorption lines. 

VK acknowledges Prof. S. Miyama for hospitality, 
and National Astronomical Observatory of Tokyo for 
providing COE fellowship. 
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